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proof of the independence of the theory under the changes of the quantum and classical gauges 
is given. The singularities appearing in the perturbative expansion are eliminated by properly 
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O . 

O ! I. INTRODUCTION 

oo ; 

(N . Improving the understanding of chromodynamics is an important issue in modern Particle Physics . Although the 
I small couplings at high momenta allow the use of perturbation theory in this region, the expansion is however unable 
^ . to describe the physics at low energies. In former works we have been considering the development of a 

Qv^ modified perturbation expansion, searching to describe at least some low-energy properties. The proposed scheme 
\^ . conserves colour and Lorentz invariances. Indeed, it was motivated by the aim of eliminating the symmetry limitations 
' of the earlier chromomagnetic field models S S ^| ■ The Feynman expansion proposed in Ref. included in 
the gluon propagator a term produced by the condensation of zero-momentum gluons. Later, in Refs. llll | the 
' interest of also including quark condensates was argued. Just in the first approximation, the modified rules produced 
J;0 [ a non- vanishing gluon condensation parameter {g^G'^) Afterwards, in Ref. 0, a non-zero tachyon mass for the 
■ gluons was evaluated at the one-loop level. This outcome was consistent with the conclusion in Ref. ^3 about the 
' tachyon character of the normal Green functions in the presence of a gluon condensate Also, the one-loop 

effective potential as a function of the condensate parameter indicated its dynamic generation. The plausibility of 
O ^ modifying the Feynman rules of QCD, including gluon condensation, was also foreseen in Ref. [l^ . In this work a 
D filling of gluon states for all momenta k < Aqcd, was assumed in defining an alternative free vacuum state, upon which 
to connect the colour interaction. 
I The addition of a delta-function term at zero momentum in the gluon propagator had also been investigated 
in Refs. [TtL ll^. The aim was to develop a phenomenological model for meson resonances, showing a confining 
inter-quark potential. The multiplicative constant of the delta function had the opposite sign with respect to the 
one chosen in Refs. 0, |^ . For this selection of the sign, the singularities of the quark propagator showed no pole on 
the real axis Thus, the results led to a physical description different from that in Refs. 0-13 ■ 

The expansion examined in Ref. was derived in Refs. 0, as the result of the adiabatic connection of the 
colour interaction, upon a particularly constructed free physical state. The operator quantization scheme of Kugo 
and Ojima |lSj.20j allowed this state to be determined by the creation of zero-momenta gluons and ghost particles. 
The propagators considered in Ref. 0, then emerged as generating the Wick expansion. The annihilation of the 
modified vacuum by the action of the Becchi-Rouet-Stora-Tyutin (BRST) charge assured its physical character at 
zero coupling. It also followed that the parameter describing the condensation should be real and positive. This 
specific definition, in our case, discarded the possibility of fixing the sign of the factor multiplying the delta function, 
in the way chosen in Refs. 1^ should be added that the propagators employed by Munczek and Nemirovsky 

were also obtained in Ref. |2ll |. as generated by the Wick expansion from a vacuum including condensed gluons in 
a squeezed state. It may be that the BRST physical condition also could be satisfied by that class of states. Then, 
their alternative approach could be constructed along the same lines. The difference in the end results for these two 
treatments might be related with the fact that two different kinds of 'squeezed' gluon states were employed [21} . 

In the present work, we consider the investigation of the gauge-invariance and regularization properties of the 
modified expansion under study. Our main objectives will be: a) to generalize the procedure to arbitrary values of 
the gauge parameter a, b) to eliminate the infrared singularities appearing as a consequence of the delta-function 
structure of the new propagators, c) to present the formal proof of the Ward-Takahashi-Slavnov identities and 
the gauge-parameter independence of the physical quantities, and finally d) to check the satisfaction of the gauge 
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independence of the effective action and the transversahty of the gluon self-energy at zero values of the mean fields, 
to second order in the coupling g and all orders in the expansion in the condensate parameters. 

Section 2 starts by presenting the initial state to be used in the Wick expansion formula. It is then employed 
in the derivation of the modified propagators. At this point, the invariant functions appearing as the results of 
the commutators defining all the propagators are modified by implementing the Nakanishi procedure for the infrared 
regularization of these functions |23| . This procedure allows us to show that all the Feynman parts of the propagators 
are regularized to be equal to zero in a small neighbourhood of the origin in momentum space. The size of the vanishing 
region is fixed by the Nakanishi regularization parameter a. This is a basic outcome that further allows the elimination 
of a large class of the singularities generated by the delta-function character of the modified propagators. 

Section 3 is devoted to deriving the Ward-Takahashi-Slavnov (WTS) identities corresponding to the Green functions 
in the modified expansion [2^0, |2^. They fully coincide with those in the usual theory (PQCD). This outcome 
is suggested after considering that the modified propagators are allowed inverse kernels of the tree-level inverse 
propagators of massless QCD. The generating functional of the Green functions is transformed into a Feynman 
integral form, only differing in the boundary conditions on the fields, from the standard one in PQCD. 

In Section 4, we formulate a procedure for the elimination of the infrared singularities appearing in any order of 
the perturbative expansion due to the S{p) form of the condensate terms. The properties allowing the cancellation of 
those divergences are two. The first is the use of dimensional regularization for the Dirac's delta function evaluated at 
zero momentum. As was shown in Ref. '26], in the case of spatial arguments, these expressions can be made to vanish 
in dimensional regularization. An identical proof is employed here to argue the same conclusion for momentum-space 
arguments. Using this rule, a large class of the singular terms are removed. The second recourse leading to a vanishing 
result for the remaining singular terms is the Nakanishi infrared regularization described in Section 2. Specifically, 
the equality to zero of the Feynman contribution to the gluon, quark and ghost propagators in the neighbourhood 
of p = (of a size fixed by the Nakanishi parameter a) allows us to cancel the remaining singularities. Since these 
singular terms arise from the evaluation of the Feynman propagators at zero momenta, this is simply obtained after 
setting the absolute values of the momenta at which the condensate states are created, being sufficiently smaller than 
a. 

Then, Section 6 presents the checking of the transversahty of the polarization tensor and the gauge-parameter 
independence of the effective action evaluated at zero value of the mean fields. This is done here up to second order in 
the coupling constant and any order in the condensate parameters. Afterwards, comments on the results and possible 
directions of their future extension are given in a summary. 

Finally a resume of the Kugo-Ojima quantization procedure for free massless QCD for general value of the gauge 
parameter a, and a complementary calculation, are given in Appendices A and B, respectively. 



II. MODIFIED FEYNMAN EXPANSION 



A. The initial vacuum for the Wick series 



Let us consider in what follows the selection of the initial vacuum state of free QCD upon which the interaction 
will be connected. It will be a generalization of the previously defined one in Refs. [^|^, but physically equivalent to 
it, since its physical part (the transverse gluons part) is not modified. The changes appear only in the non-physical 
sector, which we choose the most general possible in order to have a wider freedom in the selection of the resulting 
propagators. The main elements of the Kugo— Ojima operator quantization of the free massless QCD for an arbitrary 
gauge parameter a are reviewed in Appendix A. The resume collects the equations of motion and commutation 
relations between the fields and the creation and annihilation operators. The definitions of the invariant functions 
and the various sorts of wave packets to be employed below can also be found there. The vacuum state will be defined 
in the form 



v?) = exp <j E E 

a,=l pi^\pi\=P<a 



E 

A=l,2 



Ca {P) .a+ ,a 



At+.At+, + C3 (a. P) ( Bt+Ab''+ + ict+ct+) 

Pi, A Pi, A ^ ' ' y Pi Pi Pi Pi J 



+D (a, P) K+S^+J + E E \\^)- (1) 

s— 1,2 Pi ,qi .\pi .qi\—P <,(j 



The gluon (first part of the exponential) and quark (second part of the exponential) pairs defining the state are created 
for a set of small but non- vanishing momenta pi in order to avoid the difficulties occurring if they are directly chosen 
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to have zero momentum. A sum over the set pi for \pi\ = P < a \s introduced, also to have freedom in eliminating 
any preferential direction in the space, after the Nakanishi parameter will be taken in the limit ct — > 0. 

In order to simplify the exposition, we evaluate gluon expressions below for a generic mode of given colour a and 
momentum pi indices. This can be done because for the free theory, and for the specific quantities to be evaluated, 
contributions of different modes can be worked out independently, thanks to the commutation relations HA12II . At 
the necessary point in the discussion, all contributions will be included. 

The proof that the proposed state (Q) satisfies the required BRST physical-state conditions 



Qb\^) 



k k 



= 0, 



0, 



k,a 



proceeds as follows. 

Using the commutation relations (|A12|I . the state | ^P) in Eq. |^ can be written in the form 



*) =exp(i?(a,P)i?i+i?^+) 1$), 



where | $) has the form 



exp 



Ci (F) At+,A^, + C2 (P) At+,At+, + C, (a, P) 
and satisfies the physical-state conditions, as was probed in Refs. 0, 1^ 

Qb I *) = Qc I = 0. 
The commutation relations (jA12|) also allow to show that 



— a+ a+ 
IC- C- 

Pi Pi 



0), 



Hence 



Qb I *> 



QB,exp (^D {a,P)Bl+B'i-^ 



e^^[D{a,P)B'i+B'i+]QB 



Qcexp (d {a,P)B-+B-+) 



0. 



0, 



Qc I *> = exp (d (a, P) B^+B^+) Qc \ ^ 0. 



The state ^ is then a physical state of the free theory. The non-physical sector of | is undetectable in the physical 
world, and can neglected when calculating physical quantities like the norm, energy, and particle number. 

As stated before the physical sector of Eq. ^ is the same proposed in previous work 0,13 , then for the calculation 
of the norm a similar result is obtained 



^=(*i*> = n n n 

A=l,2 a=lp'i,|p-i|=P<CT 

The normalized physical vacuum state is defined as 



2m 



(2m)! 
(m!)^ 



for \Cx{P)\ < 1. 



I^) = ^l*). (2) 

The modifications introduced in the usual perturbation theory by the state |^ *re calculated in this section for 
arbitrary values of the gauge parameter a. For this purpose the expressions derived in Refs. |^ for the Green 
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(5) 



functions generating functional of a free theory are used 

^0 [j, 1 1 V, V] - ^o' b'] ^0 h v] zt [U], (3) 
Z'oIj] = (^-iTcxp (^i J dx r (x) (x)^ |§) = (§1 exp (^z J dx f {x) A' {x)^ exp [i J dy f {x) A+ {x)^ |§) 

X exp (^J dx dy 0{yo - xo)f (x) [A- {x) , At (y)] fiy)^ , (4) 
Zq [ry, ^] = (^1 exp J dx rj {x) il^^ (x) + ifi (x) rj (x) 

X exp (^J dx dy T]{y) 6'(a;o - yo) {V^" (y) , - 6l(yo - xo) jv^^ (y) , -0 77(2;) j , (6) 

= ZrhvK^^[rj,r]], (7) 
Zt[^,I] = {^\exp(^i J dxl^{x)c~ {x)+c- {x)^{x)]^exp (^t J dy[^{x)c+ {x)+c+ (x) ^ (x)]) |§) 

X exp (^J dx dy [e{xo - yo) {c" (y) ,c+ (x)} - 0(yo - 2:0) {c+ (y) ,5" (x)}] ^{x)^ , (8) 



exp [ i dy 



— + 



^ (x) ■0+ (x) + ■0 {x)ri{x) ]\^) 



(9) 



where the colour and spinor indices are implicit. As is clear from the above formulae, the changes in the gluon, quark 
and ghost generating functionals with respect to the standard ones in PQCD are given by the factors: 



^o''"b] = (*| exp z / f [x) A- {x) dy (x) A+ (x) 1*), 



Zq''"[77, 77] = (^tl exp ( i / dx ri{x)ip (x) + ip {x)ri{x) ) exp ( i / dy rj (x) tp^ (x) + tp~^ (x) rj (x) ]\'^), 



e,e]-(*|expU 



dx[^{x)c- {x) + c- {x)({x)]^exp(^t J dyl^{x)c+{x) + c+{x)^{x)]^ \^), (10) 



which all reduce to unity in the case of the standard PQCD. 

Let us now explicitly state the infrared regularization rules of the expansion to be considered: 

a) The invariant functions D and E, in terms of which the commutators (anti- commutators) defining the standard 
Feynman propagators are expressed, will be assumed to be given by the Nakanishi infrared regularized expressions 



E^^^{x-y\<j) ^l^—({xo~yo)jr^ r - l) \im Di^^x - y \ a,m) . 



(11) 



b ) For the evaluation of the changes in the propagators produced by the condensation, which are physically related to 
the zero-momentum modes of the theory, the unmodified field operators will be used. In addition, the spatial momenta 
for which gluon and quark condensate states are created will be assumed to lie well within a small neighbourhood of 
p = 0. The radius P of this region will be chosen much smaller than the Nakanishi parameter a in Eq. 

Let us consider in what follows the evaluation of the various Feynman propagators following the above prescriptions. 



B. Regularized Feynman propagators 

1. Gluons 



The Feynman propagator for the gluons is determined by the commutator between the creation and annihilation 
parts of their field operators. After substituting l|A3|l in this commutator, the expression can be transformed in 
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the following way 

[A^- (x) , Al+ {y)] = 5'^\g,.D+ (x -y\0)-{l- a) {dlD^^^^ {y) aliv) IklA^) + 

'k 

+d:D'--^Hx)Y,g,{x)f,lAy))), (12) 

~k 

= lim {5'^\g^..D+ (.x - y |(7, 0) - (1 - a) dld^ E+{x~y\a))), (13) 

where it has been defined D± {x — y\Q) = limo-^o D± {x — y\a,Q) . In the last line of this relation, in which the 
commutator has been expressed in terms of the invariant functions, the Nakanishi regularization pi|l has been 
introduced. The functions taken for a finite value of cr, will be employed in what follows, defining in this way the 
regularized propagators. 

Substituting the above formula in the quadratic form in the gluon sources defining the Feynman contribution 
Zq'^ for the generating functional in wc now have the form 

5t = / dxdy 9iyo - x^)f (x) [A~ (x) , A^ iy)]riy), (14) 



I I dxdyr'^x)g^Aeixo-yo)D+{x-y\a,0)+9{yo-xo) D+{y-x\a,0) )r'"(y)- 
dxdy d-^r^ixmxo - xo)E+ {x-y |a) + e{yo - yo)E+ {y - x |<T))a^j-''^(y), 



2 

(l-a) 



1 1 dxdyr^-{x)DF^,{x~y\a)r-iy)+ 
^^""^ I' dxdy dir^{x)EF^,{x~y\a)dy3^^'^{y), 



where, in order to shift the derivatives to act on the sources, the following properties have been used 

dldliEpix - y \o)) - dldliOix^ - y^)E^{x - y|a) - ^(yo - x^)E_{x - y |a)), 

= (^(xo - yo)dldlE+{x - y k) - 9{yo - x„)dyd;EAx-y\a)), 

in which we had defined 9^ — + u^9o with = {0,81,82,83), = (1,0,0,0), and the final equality follows 
because the delta- functions evaluate the factors [E+(x — y)+E^(x — y)) at xq = yo, where E+{x — y\a) — -'E-{x~y\a) 
is satisfied. Then, the Feynman propagator contribution to the generating functional of the gluons takes the form 

exp{SD + Se) = exp ( / dx dyOiyo - xo)j{x)[A^ (x) , A+ (y)]j(2/) I , 



exp Q J dxdyf^''{x)g^,^DF{ 



x-y\^)r''{y) ) X 

exp ( + / dxdy 8lf^^^{x)EF{x - y ^f^^^iy) ] , (15) 



where Df{x — y\a) indicates the Feynman term in the modified gluon propagator for a = 1, although it differs in a 
constant multiplicative factor from the usual definition. 

The generating functional expressed in momentum space for the gluons can be obtained now by passing the 
quadratic integrals on the sources to Fourier representation. The following convention for the Fourier transform of 
any quantity Q, and the Fourier transform of the Heaviside function will be used 



dQ f dQ 

Q{x) ^ I j^--^Q{p)e-Kp{-iq.x) = J j^--^TT{Q{x)\x,p) exp{-iq.x), (16) 



^(2^0)= / ^ ^—exp{~iqoXo). (17) 

J -00 27r go + «e 

Let us consider first the second term in the argument of the exponential Se, which is the most involved one. After 
substituting the Fourier representations (|15|) for the sources, the Heaviside function (|17|l and (jA2p . we can obtain 

Se = J ^fc^j^(-fc)(.Fr(0(^o)i?+(^|a)|z,fc)+.Fr(0(zo)i?+(^|fT)|z,-fc))fc^j:(fc). (18) 
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The first Fourier transform in the integrand of (|18|l . after substituting the definition of the E function and the 
Fourier representation (|17|) of 0(zq), can be obtained in the form 



TTmf W M ^1 i.^ ^0{\k\~<j) (2|fc|-fco) 

-4|fc|3(fco-|fc|+^e)2 

Since the second Fourier transform in the integrand of (|18() is simply the first one, but with the integration 
momentum k argument replaced by — fc, adding the two terms leads to 



2 J (27r)4 ^ ^4\k\^{ko-\k\+ie)^ 
-4| k \^{-ko - I k I + ie)2 

For the terms in (|18|1 associated to the invariant function , a similar but much simpler calculation thanks to 
the absence of dipole fields, allows us to obtain 

which inserted in (jSJ gives, for the regularized gluon Feynman factor in the free generating functional, the expression 



2. Ghosts and Quarks 

In a similar way, the regularized ghost and quark propagators (the quark taken in the massless limit m ~> 0) can 
be obtained in the forms 

ZI'^[V,V] = exp J j^vi-k) \j _^ J vik)j , (19) 

Therefore, the Nakanishi regularization leads to the vanishing of all the Feynman propagator terms within a tube 
containing A: = and defined by the momenta having modulus of their spatial component smaller that a, for arbitrary 
values of the temporal component fep. In this way, as will be seen further in the discussion, the adoption of this 
regularization procedure naturally leads to the cancellation of a great number of the singularities in the perturbative 
expansion. 



C. Condensation-induced changes in the propagators 



The effects on the propagators of the gluons and quarks created by the condensation of gluons, quarks and ghost 
in the free ground state, will be considered in this subsection. Let us separately examine the two cases. 
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1. Gluons 

Taking into account expressions ljA3|l . the gluon annihilation and creation parts in (|10|) can be exphcitly written as 



fc \A=1,2 ' / 

(^) = E 1 E ^tJ^:^^ + 4"^ a* (^) + [a* (^) + (i - «) (a* (fco, 2:0) z^* l,,, (x) + i?; (fc, x))] 



where A (fco, xq) and (fc, a;) are defined as 



A (fco, So) = {ikaxo + 1/2) , (fc, x) = ^^"^^° gfc(x). 

2 fc 2 fc 

On arriving to (|21|l . the term of 9^13 5^(0;) in ljA3|l was evaluated by using the wave packets (|A7|) and polarization 
vectors (|A9p with the result 

d^D[l^^^^gkix) = 9^1/2(v2)-i(xo9o - 1/2)^^^ exp (-ikx) = A (fco, xo) A.l,^ (x) + (fc, x) , (22) 
which defined A (fco, xo) and B^ (fc, x). In a similar way, for 9^Z3|^^^-' 5^(0;) we obtained 

d^nfj^^^glix) = A* (fco, xo) flL,, i'x) + B; (fc, x) . (23) 

As commented in Appendix A, the transverse term is not modified by the generalization of the theory for arbitrary 
values of a. Then, its contribution is the same as calculated previously j^, It is important to recall that thanks to 
the diagonal block structure commutation relations (|A12p the calculations are done by decomposing the expression 
(|10|l in products of contributions coming from separate modes associated to the diagonal blocks. 

The result obtained in Refs. 3, 5] for the transverse-mode contributions, for \Cx {pi)\ < 1, is 

^(0 I exp I 5] ^i^A%^,A%^, I V / ^^-^"'^ E ^ptxfp*, (^) [ 

\A— 1,2 / I A— 1,2 I 

X exp M dxJ^^'^ (x) J2 ^kxfL^^ (^) [ I E 

[ A=l,2 J \A=1,2 



CxiP) 

2 "^ft,A^p,,A 



. , -^^^.A^^A I I 0), 

= exp - ^ J« A L , (24) 

Att2 2[l-\C\iP)\') J 



where the normalization factor has been cancelled. It should be recalled that the momenta pi are very small and the 
sources are assumed to be located in a finite spatial region. The following simplified notation was also introduced in 

m 

The calculation of longitudinal and scalar-mode contributions is, on this occasion, more elaborate, and a brief 
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exposition of it can be found in Appendix A; the result obtained for IC3 (a, P)| < 1 is 



exp 



dxdy 



C3 («, P) + C3*(a,P) + 2 |C3 [a, P)\' 
l-|C3(a,P)|- 



X ( es,p^ (pi) eL,u (pi) + ^] ,^^2 l-'^L.fj^ {Pi) + i2pio(5po] <^l,u (Pi) ] + 



4 Ik I 



D* (a, P) (C3 (a, P) + lf + D {a, P) (C3* (a, P) + 1)' 



eL.Ai (Pi) e-L,!^ {Pi) 



(1 -IC3 (a, p)r 

; (0 I exp {Cl (a, P) 4'" + D* {a, P) B^Bi) exp (^3 (a, P) + D {a, P) P|+P|+) | 0), 



(25) 



where the normahzation factor (last line) is equal to one. Now inserting expressions H24|l and (|25(l in H1Q(I . performing 
some algebraic manipulations [keeping in mind the properties of the polarization vectors defined in (jA9|) ]. assuming 
Cl (P) = C2 (P), and introducing the contribution of all momenta pi {\Pi\ = P), the following expression is obtained: 



exp 



Ci(P) + Ci*(P) + 2|C i (P)r 

2fi-ici(p)r 



C3 (a, P) + C3* (a, P) + 2 |C73 (a, P)|' Ci (P) + C^* (P) + 2 |Ci (P)|' \ 



(l-\C^{a.Pf) 



(P)p) 



C3(a,P)+C3*(a,P)+2|C3(a,P)|^ (1-a) 

—2- (PifiPii^ ~ ^piod^opii, j 



l-|C3(a,P)r 



4 Ik 



D* (a, P) (Ca (a, P) + 1) V (a, P) (Q (a, P) + if 



l-\Csia,Pty 



(26) 



It is important to notice at this point that the combinations of Ci (P), C3 {a, P), and D* {a, P) in Eq. (|26|l are real. So 
that even if these parameters are complex, only their real parts contribute to the propagator. After this observation, 
they will be considered real in what follows. 

In expression H26() we take the thermodynamic limit {V — > 00), replace sums by integrals, 



2ir 



Pi,\pi\=P 



(27) 



and perform the integrations. 

As a second step, we proceed to take the limit P ^ 0, for that we consider that the parameters introduced in the 
modified vacuum state to be 

Ci(p)^l-^P, Ci>0, 
C3(a,P)~l-^^P, C3(a)>0, 



D{a,P)^D(a). 
After that, the result obtained for Eq. H2()|) in the limit P ^ is 



(28) 



{/■ 



exp <; / dxdy J"'" (x) J"'" (y) 
(1-a) 



(27r)' 3C3 {a) 



{2nfCi 3(27r)^ yCs (a) Ci 



(27r) 3[C;3(a)J 



(29) 
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In order to make explicit the Lorentz invariance, we fix D (a), such that the terms proportional to S^^qS^q cancel 
out. Then D (a) is determined by the expression 



Dia) 



1 

C'l 



2C3 (a) 



and H29|l takes the form 



exp ■ 




(30) 



In expression H30I) we notice that the term in the exponential is real and non-negative. This is because Ci > 0, 
C3 (a) > and a > (a most be non-negative for the convergence of the Gaussian integral in which it was introduced 
in the path integral formalism. It can be stressed at this point, that, if could be possible justifying the validity of the 
recursive solution as an analytical extension in the parameters, it looks possible to trace a connection with the approach 
of Refs. [rRIT8Ll2]| ). As a final step, we chose the parameter C3 (a) in such a way that the determined modification 
does not depend on the gauge parameter a, as it was selected in an earlier work 2]. Then C3 {a) = \KCi (1 + a) , 
with K a positive constant. 

Considering the above remarks and defining in the generating functional modification an alternative non-negative 

constant Cg = (l + -^) (which we call condensate parameter) for further convenience, expression (|5n|) takes the 

form 



exp I / dxdy V J'^-'^ {x) {y) ^ 

[•^ a=l (27^ 



From the above construction, the parameter D (a) = D (1 + a) , where D is an arbitrary real constant (equal to 
zero for K = I), and the vacuum in the zero-momentum limit has the form 



= exp 



E Uti^oS + Uo^A^,% + B^,+A^'^+ + tc^+4+ + D{1 + af BS+B^,+ 



0), 



(31) 



from which the state selected in the previous works |^, corresponds to the particular case a — 1, C3 {a — 1} — Ci 
and D (a = 1) = 0. 

Thus, it follows that the vacuum state introduced in l|31() modifies the usual perturbation theory only through a 
change in the gluon propagator which now takes the form 



DfA^ -y) = 



dk 



fc2 + 



ffAii/ - (1 - a) 



ky 



iCgS {k)g^^ 



exp {—ik {x — y)} . 



(32) 



in which the first term is the usual Feynman propagator, but now including an infrared regularization. 



2. Ghosts 

The change in the ghost generating functional is the same as calculated in Refs. This is because the ghost 

sector in the vacuum state remains unchanged. As was mentioned in these references for the value C3 (a, 0) = 1 as it 
was fixed here, there is no modification of the ghost sector and its generating functional and propagator remain the 
same as for the usual PQCD. 
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3. Quarks 

Let us now consider the modification in the quark propagator produced by the presence of its condensate. The 
factor corresponding to quarks in the generating functional can be written as 



^o''"[^,^] = (0|exp 



J2 J2 4.^1-. I exp<{i / rj{x)yj {x)+i> {x)r]{x) 



X exp <i dy 7] (x) ip'^ (x) + ip (x) rj (x) 



s=1.2pi,gi,|pi,gi| = P<(T 



= (0 I exp [Cfj:,a''fbf,] exp \^~crfaj + ayfe+j exp {-cr}a/ + a*fbf} exp |^C//'feyta^+ 




(33) 



= ^o'" K ^] = ^o'" M , V = (ct, a), I (?) = exp 



|0), 



where the quark-field operators were expressed as sums over the creation and annihilation components. In order 
to simplify the discussion the following compact notation has been employed /, /' = Pi,i — 1,2,3, N. The new 
indices, / and /', correspond to the values of the momenta for which condensate states are created. The spinor index 
s = 1 or 2 have been omitted, since the discussion can be done almost up to the end for each value of s, thanks to the 
commutativity of the creation operators for different indices. The coefficients of the quark creation and annihilation 
operators and b"^ are expressed as 



dx u~ ri{x) exp(ifca;). 



dxri{x)v^ exp{ikx). 



(34) 



where u and v are the usual spinor solutions of the Dirac equation defined in Appendix A. 

Let us now define the transformations 11^^ = exp(— a/a^ + o'/6y^) and U2 = exp(— cr^a/ +W*j:bf), which acting upon 
the creation and annihilation operators through a similarity, lead to 

UibfUjP^=bf -Wf, UiafU^^ = aj + af. 



U2b+U2^- 



+ (7 J, UiOjUi ^ = Oj^ — a J. 



After also consider {0\Ui = (0| and C/2IO) = |0); it follows that 

Z„''™[v] = (0 I exp [C} f,{af + af){bf, - a^,)] exp [Cf f,{b'+ +a},){a'+ - a})] \ 0), 

= {q\ exp [-C} f,(Tfaj, - Cff>a*f,a}] exp ^-a/ Cff>af> + afC^j,bf,\^ exp 

It is now possible to employ the following identity for operators that commute with their commutators 

exp(A) exp(B) = exp(B) exp(A) exp([A, B]), 



(35) 



-a)Cjf,a),+<j*C^^,b+ 



(36) 



to change the order of the exponential operators appearing in the mean value on the quark 'squeezed state in Eq. 
(|35|l . The commutator between the two arguments of those exponentials can be calculated to be 

-^afCffiOf, +afCjj,bf,,-a*fCff:a), + a}C ^ ^,b+\^ = - [Wf>Cj*fCff"Wf,, + ajC} f, Cj, j„a}„] , 

which, after being considered together with (|36|l . allows us to write the quark modification as 



exp 
exp 



-crCa - CT*Ccr* - aC'^Ca* - aC^C 



(37) 



1 



1 



^viRvi - -v*i?vj + viR V* 



Zq[RICvi]. 



where it was assumed that Cff is real, and the following matrix notation was introduced 
aCa^afCff,af,, {C^)ff, = {C)f>f, vi = d | , R 



-C 
C 



11 



and /C is simply the complex conjugate of a vector expressed in matrix notation. Relation (|37|l . after being applied n 
times, leads to the following expression: 

Zr[v] = expjivii? (j^ R'A VI - UlR [j2 R'A +vii?2 (j2 R'A vA Z, [(/Ci?)"+%i 

I. \m=0 / \r?j.=0 / \m=0 / J 

However, the sums of powers of the matrix R can be written in the form 
which gives in the limit n — > oo, assuming that all the eigenvalues of the symmetric matrix R^ are smaller than 1: 
Hence, the modification of the quark propagator takes the form 



where the r.h.s. of the generating functional is zero because lim„^o(^"^^) = 0, if all the eigenvalues of R^ are 
assumed to be smaller than the unity. For the argument of the exponential above it follows 

i?2 ^ dd^ , _ c^c 

vi-; ;^v, — (T a —a — —a . 

Now, assuming that the matrix C is antisymmetric and also that it satisfies C C = KI (with K including also 
possible negative values), it follows that the sum of the first two terms in the argument of the exponential in H38|) 
vanishes. The one remaining take the form 

R * ^ ^ * , — *\ ^ I * I *^ ^ ST^ ( s s* I s— s* 

1 - i?2 1 1~K^ ' i-K^ ^ f ^ i-K ^ \ I' k ^ k k 

k 

where the index s — 1,2 takes the fixed value that has been used for the whole evaluation. Now employing the 
definitions H34|l of a and a and the identities 

s=l,2 s=l,2 

and assuming that the momenta arguments pi are taken at their vanishing value in the limit L — > oo, it follows that 

K m f _ 

= - ^ _ j^ Y^ / dx dyq{x) ri{y) e^p{ik{x - y)), 

where the number N of the momenta for the created particles and antiparticles appears in the last line from the 
assumption that for small momenta the integral appearing is almost independent of fc. Fixing K as given by 
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and considering (|19|l allows us to write for the change in the quark generating functional and propagators 
Zo"'[n,v]=exp (^Cg J dxdyf]{x)'q{y)^ Zq[0], 
Zo^[v,v]=cxp (^J dxdy l]{x)G f {x - y)'q{y)^ Z,;[0], 
Z«[77,77] =exp { f dxdy rj{x)G''{x - yMy)) Zg[0], 
ie{\k\-a)-f''K 



G'^{x-y) = 



dk 



fc2 + ie 



+ CqI5{k) 



exp{—ik {x — y)). 



(39) 
(40) 
(41) 

(42) 



III. REGULARIZATION OF THE SINGULAR TERMS 

A main technical difficulty for the implementation of the expansion proposed in Refs. |^ ^ IH is that 
the modifications of the usual free propagator terms (to be called below the 'condensate propagators) are given by 
Dirac's delta functions of the momenta. This circumstance produces singular Feynman diagrams in the perturbation 
series even after dimensional regularization is introduced. The singularities correspond with the appearance of delta 
functions or standard Feynman propagators evaluated at zero momentum, after some loop integrals are performed. 
These factors occurs because of the momentum conservation in the vertices of the diagrams. Let us consider a vertex 
having n legs (n = 3, 4 for QCD). Then, when n — 1 different condensate lines join to it, the momentum conservation 
forces the value of the momentum at the only remaining line to vanish. Therefore, if a condensate line is attached 
to this ending, delta functions evaluated at zero momentum will appear. On the other hand, when a usual Feynman 
propagator is connected to this last leg, a factor equal to its value at zero momentum appears. This situation should 
be solved before a sense could be given to the modified expansion. Below, we propose two ways for the elimination of 
these singularities, which clearly should be the subject of further examination about its consistency. Let us separately 
consider in what follows the two cases. It will be assumed that the Feynman diagrams are constructed in D dimensions 
and for a fixed value of the Feynman parameter e. Before advancing, let us precisely set the rules for the order of the 
limits associated to the various regularizations that have been done up to now: 

a) The large volume limit will he taken by taking the size of the quantization box L going to infinity. This limit will 
lead to continuous values of the momenta. 

h) The Nakanishi infrared regularization parameter a will be chosen as a function cr(_L) of the spatial size of the 
system L, vanishing as L — > oo. This constant sets to zero all the Feynman propagators for the spatial momenta lying 
within a sphere of radius terms cr, and all values of the temporal component fco- 

c) The maximal size P of the set of the spatial momenta ki for which condensate states are created will he assumed 
to he another function P{L) also vanishing in the limit L oo. However, it will he chosen smaller than cr{L), in a 
proportion to he defined in the following. 

d) Finally, the auxiliary mass m of the quark field will be chosen as a function m[L) tending to zero for L oo 
with a behavior to be also described below. 

The limit L — > oo, in which <t{L), P{L), and m{L) vanish, will be taken in the first place. With this step we recover 
the Lorentz invariance of the Feynman diagrams that continue to be functions of the dimensional regularization 
parameter D, it mass scale /i, and the Feynman constant e. 



A. S{0) singularities 

An idea that comes directly to the mind when considering the singular terms of the form 6(0) is that they should be 
considered in dimensional regularization. But, it has been argued that the delta- functions evaluated a zero spatial 
coordinates can be analytically extended to continuous dimensions D and, moreover, that their expressions vanish 
in the limit Z) — > 4 |2a |. This is not necessarily an un-natural result, as it will illustrated in a second method 
to remove these singularities, already in the step when the limit L — s- is taken. Helpful in this purpose will be 
a proper selection of the set of momenta {ki} for which the condensate states are created. Let us start, however, 
with the dimensional regularization argument. We follow the same procedure here and interpret the Delta functions 
appearing as D-dimensional ones. The Wick rotation of the momentum integrals will be assumed to be already done. 
Then, it is possible to reproduce, step by step, the arguments of Capper and Leibbrandt to conclude that these 
factors should vanish after removing the dimensional regularization. Let us do it below for the sake of concreteness. 
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For S{Q) we can write 



S{0) 



This is a singular Z3-dimensional integral in Euclidean momentum space, similar to those in real space, considered in 
[2^. Then, it can also be written as follows 



dp 



D 



{2n) 



dp^ p^ 
(2^ 



dp^ 

(27r)^ 



exp(— 



However, using the redefinition of the generalized Gaussian integral for continuous values of the dimension D con- 
structed in Ref. (26j . it is possible to write 



dp 



D 



dpD ^ 
IT 



and 



(2^) 



m 



exp(— sp^) : 
sp')- 



1 



{4TT) 



■ exp 



D 



exp(. 



(47r)- 



ds- 



{Any 



D 



f 



exp 



exp 



-sf 



where / is the function introduced in '2S\ for extending the generalized Gaussian integral formula to non-integral 
dimension arguments. These functions vanish for all integral values of D. Then, after using the integral definition of 
the Gamma function r(2;) = dt t^~^ exp(— t), it follows that 



6{0) 



(47r)"? 



D 



Til- 



which vanishes exactly in the limit D ^ A. Therefore, we will assume that the evaluations associated to the modified 
expansion will be made by using the above representation for the factors 5(Q). As a consequence, it will be considered 
that all the diagrams in which such a kind of singularities appear will vanish in dimensional regularization. 

Let us consider now a property through which it seems that we can directly eliminate such singularities, even in the 
absence of the above exposed dimensional regularization argument. The main observation is that all the 4- momenta 
of the modes associated to condensate states are null 4- vectors lying in the forward light cone. Therefore, if we define 
the set {k i\ of the spatial momenta components used in the construction of the 'squeezed states, with the unique 
condition that all of them have different orientations in their spatial parts, the sum of any subset of three or four 
4- momenta can not add up to a zero result, since the temporal components are by construction all positive in the 
forward light cone. Therefore, the imposition of the rigorous conservation of the 4-momentum can not be satisfied 
by the sum of any set of three or four momenta arriving to a given three or four legs vertex. Hence, this observation 
indicates that the elimination of the 5(0) singularities is in fact a natural result of the way in which this massless 
problem has been regularized, which leads to the same outcome given by dimensional regularization. 



B. b) Feynman propagator at p = singularities 

For the elimination of this kind of singular behaviour, the Nakanishi infrared regularization implemented on the 
Feynman propagators in a previous section is specially helpful. These singular terms are simply vanishing because of 
the appearance of the infrared regularizing factor 0(|"p | — o) in each of the Feynman propagators, before taking the 
limit L — > 00. 

Let us examine below the loop expansion after assuming the above prescriptions and the modified free propagators 
for a gauge parameter a : 

DtM m) = '-^-^ - _ S^^Sip), (43) 

"'^ p'' + le \ p'' + ie J 

^ p^ + le 

'~'gh\P) — J 
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and the standard vertices of QCD. Then, it follows that all the diagrams having a fixed number of loops showing 
both types of singularities will vanish in the limit L —f oo. Therefore, after taking this limit, the remaining finite 
diagrams (due to the dimensional regularization) can be evaluated by using non-distorted dimensionally regularized 
propagators, following from l|43|l in the limit cr — > 0. 

Assuming that the chosen conditions have eliminated all the singular contributions, it is possible to comment about 
some properties of the diagram series based on the above propagators. For example, it follows that the appearance 
of a number m of the condensate propagators within an n-loop one-particle-irreducible(lPI) diagram will eliminate 
m of the n-loop integrals associated to this contribution. Therefore, the considered diagram will now be an 'effective' 
(n — to)-1oop one. Consider also, expressing the condensate parameters Cg and Cq in favour of the new ones: 



which incorporate a power of 2 factor of the coupling constant g. Then, the 7i-loop IPI diagrams of the effective 
expansion, considered as power series in the new three parameters m^, Sf and g, also show the property that, given the 
number of external legs of the diagram, the number of loops is fixed by the power of appearing in it. This conclusion 
directly follows from the fact that each time a condensate internal line appears, a loop integral is annihilated, and 
correspondingly the power of g of the diagram is reduced by 2 in the new expansion. Let us consider that Pg is the 
power of g corresponding to an n— loop IPI diagram in which the parameters are the original ones. Thus the new 
power of g of this diagram with m condensate lines when the new parameters are introduced will be 

= - 2m. (45) 

It seems that this property, can allow for a useful reordering of the perturbation expansion. To see an indication for 
this, let us consider a particular n-loop IPI diagram in which the change of the parameters (|44|l have been introduced 
and corresponding line symbols have been assigned for the standard and the condensate propagators separately. Then, 
for any particular standard type line in this diagram, consider the infinite summation of all the zero order in g (tree) 
contributions to the connected propagator. Any of the added terms, by construction have the same number of loops, 
but comes from higher loops in the original expansion. This is done by considering fixed the other standard lines. 
Thus, if not stopped by any difficulty associated to the combinatorial and symmetry factors in the diagrams, this 
infinite addition seems can be done for all the normal lines. The resulting modified (zero order in the new expansion) 
propagators are no other things that the ones employed in Ref. i^i] in the gauge a = 1. Therefore, it seems possible to 
demonstrate that the loop expansion can be reordered to produce an alternative version with modified propagators. 
The investigation of this possibility will considered in next works. 



IV. WTS IDENTITIES 



This section is devoted to showing that the generating functional constructed with the modified propagators, 
including condensate terms and the gauge parameter a, satisfies exactly the same WTS identities as those associated 
to the usual PQCD. Let us consider for the cited purpose the complete generating functional of the theory, which 
can be written in the form 

( / q9 s:3 qg <:4 

r' Sr]i6JkS{-i]i) v' S^iSjkS{-^i) 

in which the concrete form of the action 5* assumed and its conventions are defined in the Appendix A, and DeWitt 
compact notation for the integration and summation over the corresponding spatial, colour or spinor indices was 
employed. The specific set of indices associated to the Latin letter j, k, I is specified in a natural way by the particular 
type of field to which it is associated as a subindex. For example in ji the meaning is i = {x, fi, a). The free generating 
functional l|2Jl can also be expressed as a mean value in the initial state |^) of the evolution operator U of the 



Zo [J, V, V,^,^, 
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interaction Lagrangian of the external sources, as follows: 



T exp i d X 



Zo[0, 0,0,0] exp 



+ ^ C +c£^ + T] 1p + TpT] 



Consider now two complete bases of eigenstates \q^) of the field operators taken at the fixed times T and —T 
Thus, for these states 



4° 



?«(±r, -?) |q±) ^ <± (^) S'(±T, 1?) - -"^ 



V;'^(±T,^)|g±) 



(±T,l?)|g^ 



The explicit construction of the field operators for gauge theory in the holomorphic representation is done in Ref. [231 ■ 
This representation is the most convenient for the construction of the basis states ,3Qj . Inserting the completeness 
relation at the instants —T and T leads to 



^0 [j, '7, ^ = lim X ^^l*^^) ~T\J, r,, fj, ^, 



(47) 



-ic»i.c''a''c + V'(«7''9,. - to)?/; + jA + <^ c + c£,+J] t/j + ipr])]{q \^), 

where 2?<i> means the path integral differential VIA, c,c,'ijj,'ijj}, and the usual representation for the matrix element 
{q'^\U{+T, —T\J, Ti,f],£^,^)\q~) of the evolution operator between eigenstates of the fields has been used |22j; Af is the 
usual source— independent normalization constant fixing Zq[0] = 1. The {q~ ,q'^) dependence of the functional integral 
is contained in the boundary conditions of the field integration variables, which should be fixed to the eigenvalues 
^g±,M ' ' '"^9* ' "^9* times ±T, respectively. 

By inserting (|47|l in H46|l . and acting with the interaction part of the Action (expressed in terms of the functional 
derivatives) on Zq, it follows that 



— exp [i ( 

SjiSjkSji 



kl 



qgh 
'^ikl 



)]x 



0" + 



(48) 
(49) 



It can be noted that performing the derivatives over the sources of the interaction Lagrangian inside the integrand 
of the Gaussian functional integral is equivalent to adopting a perturbative definition of the functional integral for 
the interacting theory, as discussed in Ref. ■ 

That is, the classical action defining the complete generating functional as a functional integral is exactly the 
same as for usual PQCD. The changes in the diagrammatic expansion are thus only associated to the boundary 
conditions at —T and +T. It should be noted that for the case of PQCD the wave functions $(?*) — (9^|0) 
are specific Gaussian functions of the field components. It can be shown that the {q^\'^) are also Gaussian wave 
functions, but defined by a different quadratic form as the argument of the exponentials. 

Therefore, it can be concluded that all Ward identities of the usual PQCD should also be satisfied in the modified 
expansion. Since the action is identical, it therefore has the same BRST invariance transformation: 



sc-{x) = -Isxr'^c^xyix)., 

S^pix) = ig6Xc''{x)Taip{x), 



Sc^ix) = -iSXa^Af'-'^ix), 
5^{x) = -ig5\c''{xy4^{x)Ta, 
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which leads, thanks to its Jacobian, equal to the unit, and after applied as a change of variables in H49|l . leads to the 
functional form of the WTS identities: 




T:^il\x)-igr^^{x) 



-T, 



i5£,''{x) " iSlfix) 



iSjjlix) 
^Z[j,rj,fj,^,^ =0. 



The satisfaction of these relations, implies that the mean values for physical quantities should be gauge-invariant, as 
was demonstrated in general form in Refs. 0, . 

It should be remarked that the proof of the gauge invariance of the S matrix heavily relies on the existence of a 
mass shell defining the asymptotic states for the particles. But, precisely, QCD is expected to show the absence of 
these exact mass shells for its elementary fields: the quarks and gluons. As, has already been evaluated in Ref. Q, 
these mass shells for quark and gluons are not following in the first approximation in the condensate parameters and 
coupling constant. 

In ending this section, it can be noticed that the above outcome was suggested by the fact that the new propagators 
also are inverse kernels of the equations of motion of the free fields in massless QCD, as indicated by the relations 



1 



9u 



(1 - a) k'^K 



+ ie) 



iCfI5{k) 



- ig^Cg 6{k) 
1 

fc2 + 



(50) 
(51) 



V. GAUGE INVARIANCE IN ORDER g'^ 



Let us now study the validity of the general gauge— invariance properties that were obtained by means of the 
formal path— integral procedures. The new expansion now has the three parameters: the coupling constant g and 
the two new independent and dimensional ones associated to the gluon and fermion condensates. Therefore, 
the invariance properties for a given quantity should be obeyed in each order of a triple Taylor expansion in those 
constants. The validity of the a independence will be examined for the expansion of some quantities up to the 
second order in g and any order in the other two constants. In the first place the transversality of the gluon 
self-energy, which is one of the basic WTS identities, will be shown to be satisfied by the gluon self-energy in the 
above-mentioned orders. Afterwards, in the same approximation, the dimensionally regularized effective action will 
also be shown to be gauge-parameter-independent. In what follows, the wavy gluon and straight quark lines appearing 
in the Feynman diagrams without any addition will mean the complete free propagators, including the Feynman as 
well as the condensate components. The same lines including a central empty circle will indicate the Feynman part 
of the propagators (incorporating its a dependence); and the lines showing a central black dot will represent the 
condensate parts. A wavy line including a central transversal cutting segment will mean the usual gluon propagator 
in the Feynman gauge, having a g^i, Lorentz-tensor structure. Finally, a wavy line, but including a central coUinear 
segment, will represent the 'longitudinal' part of the gluon propagator, showing the Lorentz tensor structure (1 — Q;)p^ 
Pu- This term includes all the a dependence of the free propagators. 



A. Transversality of 11^,^ in order 

Figure n shows the diagrams contributing to the polarization operator up to the second order in g. Therefore, 
collecting the terms in these expressions having certain definite order in each of the two condensate parameters will 
define the expansion coefficients of the polarization tensor in the series of the three parameters. All the terms of 
non-vanishing order in the condensate parameters are associated to the last four diagrams in Fig. ^ 

The first four diagrams in this figure represent the usual second-order contribution to 11^1^ , which is known to be 
transversal. The condensate-dependent term can also be represented as in Fig. [3 

Note the absence of terms coming form the quark condensate. This is a consequence of the vanishing of these terms, 
which is directly due to the fact that a trace of an odd number of 7 matrices appears in their analytic expressions. 
Further the last two diagrams in Fig. |21were the ones evaluated in Ref. 3 in the Feynman gauge and whose result 
is transversal. Finally, the analytic expression of the first diagram can be written in the form 
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> 






FIG. 1: The diagrams contributing to the polarization tensor up to the second order in g and all orders in the condensate 
parameters 





FIG. 2: Condensate dependent contributions to the polarization tensor. 



(1 Oi)Cq I ^2-^ faia2a fa\a2b 



Till/ — 



which, after employing the function defining the 3-gluon vertex j27l | 

^AtiA'2M3(fci>fe7^3) = (fci - ^2)^133^1^2 + (^2 - h)^i^9^i3^i2 + -~ ki)fi^g^^^^, 

leads to 



T, 



ab _ (1 '^)Cq9^ raia^a raia^b ( PfJ.Pi' 



showing the transversality of the polarization tensor up to second order in g and all orders in the gluon and quark 
condensate parameters. The a independence of the effective action in the same approximation will be studied 
below. 



B. a— independence of the effective action in order g 

Figure |31 shows the terms of order 2 in the coupling constant (and any order in the other two parameters) of 
the effective action, evaluated at the zero value of the external fields. Expanding the gluon propagator as a sum of 
its Feynman gauge component plus the longitudinal and condensate parts; and the quark one as the usual Feynman 
propagator plus the quark— condensate part, the effective action can be expressed as shown in Fig. O 

Observing this figure, it can be noted that the whole a-dependent contribution is given by the first four diagrams. 
However, the first two of them vanish because the longitudinal propagator is contracted with the polarization tensor 
in the g^ approximation, which is transverse. Further, the third term also vanishes thanks to the fact that the 
longitudinal propagator is contracted with the particular self-energy part, which was shown, in the last subsection 
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FIG. 3; Diagrams representing the contributions to the effective action in the second order in g and all order in the gluon and 
quark condensates 

to be transverse by itself. Finally, the self-energy quark term in the fourth diagram identically vanishes, because of 
a trace of an odd number of 7 matrices. This completes the proof of the gauge invariance of the effective action in 
second order in the coupling constant and any order of the condensate parameters. 

In concluding this section it can be remarked that since massless QCD has no initial dimensional parameter, the 
removal of the dimensional regularization in the diagrams of the theory can be done only after partial summations 
are performed. These summations, then, will input the new dimensional parameters associated to the condensates, 
implementing in this way the dimensional transmutation effect 32] . This and other issues are expected to be addressed 
in next works. 



VI. SUMMARY 



The gauge-invariance and regularization properties of the perturbative expansion for QCD considered in Refs. 
0, 0, was investigated. If follows that the singularities produced by the Dirac delta- function form of the new 
terms added to the gluon and quark propagators, can be properly eliminated by combining dimensional regularization 
[26| with the infrared regularization procedure for the operator quantization of gauge theories 22] . The dimensional 
extension allows to cancel the singular diagrams in which factors (5(0) appears due to the momentum conservation laws 
at vertices, when all lines arriving are associated to condensate (5-like contributions. Further, the Nakanishi infrared 
procedure allows to get rid of the singularities in the form of a Feynman propagator evaluated at zero momentum. 

Those factors appear when n—1 lines joined to a vertex of n legs are of the condensate kind, and the remaining one 
is a Feynman propagator. In connection with the gauge-invariance properties, it is argued that the modified expansion 
should satisfy the same WTS identities as the usual PQCD. In addition, it follows that the formal functional integral 
representation of the generating functional of the Green functions only differs from the one associated to PQCD in the 
boundary conditions for the fields at t = ±00. The formal functional integral results are checked in the second order 
in the coupling constant and any order in the condensate parameters. Firstly, the transversality of the gluon self- 
energy in the mentioned orders is shown. Afterwards, the contributions to the effective action evaluated at vanishing 
values of the mean fields, is shown to be gauge-parameter-independent. The work is expected to be extended in 
various directions. One of them is to implement the partial summations of the diagrams, the possibility of what 
was advanced in the text. This is needed in order to make the dimensional transmutation effective, thus, allowing 
to perform calculation depending on the new dimensional condensate parameters. Another activity to be considered 
will be to specify the renormalization procedure, which in the present problem involves three parameters. Finally, 
a task of direct physical interest will be to incorporate the knowledge about the invariance properties gained in the 
present work. It could bring a better understanding the results for the effective potential obtained in Ref. [s^, in the 
Feynman gauge a = 1. The effective potential calculated in this work gave signals of a strong instability of massless 
QCD, under the generation of a quark condensate. But, it was precisely the possibility of this outcome, that was the 
main motivation for introducing the quark condensate in the initial free vacuum in Therefore, the check of the 
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re-appearance of the instability result within a gauge invariant calculation will give further support to the dynamical 
generation of quark masses and condensates within massless QCD. 
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APPENDIX A: KUGO-OJIMA QUANTIZATION PROCEDURE FOR QCD 

The main elements of the Kugo— Ojima operator quantization of the free massless QCD for an arbitrary gauge 
parameter a will be reviewed in this appendix. The formulae will be used in the construction of the modified 
perturbative expansion done in Section 2. The classical action for the interacting fields in the Kugo-Ojima analysis 
has the form 



s 


= y {£g + Cgh +Cb+ Cg) , 






Cgh 


= -id^'c''{x)Df{x)c\x), 




= '^{x){i-i^'Dl'{x) - m5''')xi;'{x 








= 5'^'d^-gr'^Al, 




= 5^^d,~gTrAl, 







The equations of motion for the giuon, quarks, ghosts and auxiliary B fields for a gauge parameter a, in the Kugo 
and Ojima quantization scheme for free massless QCD take the form 0, 

d^Al (x) - (1 - a) 9^ (x) = 0, 
d^A"^ {x) + aB" (x) = 0, 
d^B" {x) = c'^ (x) = c" (x) = {i 7^9^ - m)V'(x) = 0, (Al) 

where, for the moment, quark fields are also considered as having an auxiliary small mass m. The notation for the 
Lorentz indices and field quantities will follow the ones used in Ref. |23|. In the case of a general value of a, the 
solution for the gluon field operator is the only one that differs from its counterpart in the Feynman gauge a — 1, 
which was considered in Ref. • As derived in 0, |23| , the non- vanishing commutation relations among the fields 
have the forms 



K ix),B'{y)] 
{c" ix),c''iy)} 

D [z I m) 

s{x) 



-S-" {g^^D {x - y|0) - (1 - a) d^d^E {x - y)) , 
5'^'>{-d^D{x-y\Q)), 

i S'^'D {x~y\0), {ij (x) ,^ (y)} = 5^' {i^d^ + m)D (x - y | 0) , 

^-r- / dk sfka) Sik^ — rn?) exp(~ik.z), E (z) — —-^ 1 -^ot^ 1 | Inn D (z I m) , 

Try J 2Vj \ dzQ J m^o 



(2- 



{x)-e{~x). 



20 



The invariant function E appear in the commutation relations for the gluon field when the gauge parameter differs 
from the Feynman gauge value a = 1, since not all the gluon fields satisfy the D'Alembert equation. More generally, 
functions are associated to each of the invariant functions D, D^, Z3_ and Dp according to 



d 11 
(z) = - lim TT^Df.) {z \ m) = 



zo 



d 

dzo 



1 I lim (z I m 



where (•) means any of the subindexes of the functions D, D^, D_ or Dp. 

The gluon and Nakanishi B-field operators solving the above equations of motion are given as 0, |23| 



(A2) 



(^) = E E ^£ (^) + 4''fk,L,, (x) + Bl [h,s,, (x) + (1 - a)d,D\l/^'^gk{x) 
B-{x) = ^l9k{x) + h.c, (x) = J2 C£ 9k{x) + h.c, r (x) = c| gk{x) + h.c, 
^^^) = E exp(-i kx) + 61+ exp(+i kx)), 

k.s 

i^{x) = E Tfy '^^S'^^ exp(-i kx) + ai+ ut exp{+i kx)), 

k.s 



h.c. 



(A3) 
(A4) 

(A5) 

(A6) 



where h.c. means the Hermitian conjugate of the previous term. The wave packets g and /, the polarization vectors 

(^)' ^i.M (^)' es, 
defined as jj, 20]: 



ej^ (fc), £L,fj. es,ft (fc) , the dirac spinor u and w, and the intcgro-diffcrential operator appearing in ljA3|l are 



5fe (x) = 



1 



V^Vko 
kx — koXQ — k.x, ko 



exp {-i kx) , /f (a;) 



1 



VWko 



eZ^ (k) exp {—ikx) 



k ■ (fc) = 0, (fc) = 0, (A:) • (fe) = (5'"^, for ct, r = 1, 2, 

T -i(|fc|, fc 



2|fc|^ 



it"(0) 



The operator iD.^^(^-' works as an "inverse" of the D'Alcmbcrtian for simple pole functions 12^, that is 



a ■ k 



v^ik) = a/tti + fco ™+'^o 



w"(0) 



w''(0) 



(A7) 
(A8) 

(A9) 

(AlO) 

(All) 



d'D['J^'^f{x)^f{x) if d^f{x) = Q. 

A large cubic box of volume y = in a particular Lorentz reference frame is assumed for the imposition of 
periodic boundary conditions on the fields. Accordingly, the spatial momenta in the above sums take the values 
k = ^(711,712,713), for arbitrary integers 711,7^2 and 713. The four- vectors for all the particles but the quarks are 

null ones, k — {\ k\,k), and for the quarks k = {\J \k p + ttt,^, k). The table below shows the commutation relations 
between the creation and annihilation operators for the various fields 
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At 

fc',A' 


A^' 

k' 


k' 




—a+ 

c" 

k' 


r' + 


1 r -\- 

h' 






















. T n 

k 








-S'-^'kk' 














Bi 

k 


































































-\ 

















S"'Srk' 





'I 























> . 



(A12) 



APPENDIX B: LONGITUDINAL AND SCALAR MODES CONTRIBUTION 



The contribution of longitudinal and scalar modes is determined by the expression 
(0 1 exp (a, P) BiA^f + D* (a, P) B^Bi x 

X exp I* 1 dxJ>^^^ [x) (4'"+/* (x) + [/* ,5,,, [x) + {l-a) {A* (po, x^) /* ^^^^^ (a;) + fe, x))] 

X exp |i y" dxJ^'^'' {x) (a;) + B% [f^^^s,^ {x) + (1 - a) (A (po, xq) /p.,L,^ (x) + B^ {p,,x))] 

X exp [Cs (a, P) 5^+4;"+ + («, P) P,-+P;-+] I 0). 
As in Refs. 0,01 j we first calculate 

exp |^ j dxJ^^'' (x) (^A^ffp„L,^. (x) + Bt (x) + {1 - a) {A {po, x„) fp^^L.^ (x) + Bp {p„ x))] 

X exp [C3 (a, P) P|+4;'^+ + D (a, P) i?2+P;^+] | 0), 
= exp jca (a, P) (^Bt+ - i j dxJ^^ '^ {x) fp^,L,^ {x)^ x 

X (^r^-* / (^) [/p.,s,M (a;) + (1 - a) {A{po,xo)fp,,L,A^) + B^.{p^,x))]j 

+D{a,V) {Bt+~i J dxJ''^''{x)fp,,L,^.{x)^ j dSxJ'^-'^ (x) (x)^ | | 0). 



(Bl) 



(B2) 



A similar expression is obtained acting to the left in Eq. I|B1|I . Substituting the above results in IjBip . and introducing 
the following notation, 

C* = C3*(a,P), C = C3(a,P), D* = D*{a,P), D = D{a,P), 
A+ = A^:''+ , i = , B+ = Bt+, B = B% , 



ai = -i J dxJ^"'" (x) [/p* ,s^p (x) + (1 - a) {A* (po, xo) /*, ^i,,^ {x) + B*^ {p^,x))] , 
02 = -* j dxJ^ °- (x) [fp,,s,tj. (x) + (1 - a) (A (po,xo) /p^.L,;. (x) + (Pi,x))] , 

62=-* f dxJ^ ''ix)fp„L,p{x), 



(B3) 
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one obtains 



(0 I exp {C* (i + ai) (s + fei) + Z?* (s + 61) (b + &i) } 

X exp {C (b+ + 62) [a+ +a2)+D (b+ + 62) (b+ + 62) } | 0). (B4) 

Following exactly the same procedure previously described in Refs. 0,13, a recurrence relation is obtained for Eq. 
B4l. 



I 

exp <^ C*aibi + D*hl + [C {C*ai - 02) + 2CD*bi] {C*bi - 62) J2 ['^1 



2(2m) 



\c\ 



2(2m+l) 



m=0 



m— 



■ (C*&i - 62)' [(2™ + 1) + 2 (m + 1) + 

■ (0 I exp^^C*AB + D*BB^ 



|2(2m) 



+C^D* [^^ + (2m + 1) |C|' 

X exp {c*"+iC"+i (C*fei " 62) i + [C*"+iC"+i (C*ai - 02) + 2C*"+iC"+ii:'*fei + 
+ 2C*"C" (n + 1) {C*D + CD*) {C*bi - 62)] ^} exp |cB+i+ + DB+B+^ \ 0). 
We then assume that |C| < 1 so that in the limit n — > 00 
lim |C|^" = lim n|C|^" = 0, 



n 

hm y nci'('") + ici 



2(2m+l) 



m— 



i-\c\ 

hm y bm+l)|C|'('")+2(m+l)|C|'('™+^' 

n— i-oo ^ L 



(i-icr 



which allows to rewrite Eq. IBSjl as 



exp < 



C*aibi + D*bi + [C {C*ai - 02) + 2CD*bi 



(l-|C|^) 

(0 I exp |c*iB + D*Bi3^ exp ^^CI3+A+ + L>B+B+| | 0). 



+ (C*6i - 62)' 



(D + C^D* 

i-icr 



(B5) 



(B6) 



(B7) 



Replacing the compact notation l|B3p in (|B7|) . and expanding all functions of pi in the vicinity of pi ~ (keeping 
in mind that the sources are located in a space finite region), it's sufficient to consider only the first term in all 
expansions. After that Eq. (|B7|) takes the form (|25|) . 
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